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ABSTRACT 

The purpose of this thesis was to postulate a method of long-range 
inertial navigation, using only acceleration inputs, and to examine theo- 
retically the response of this system to external disturbances, using vari- 
ous types of mechanizations. The problem was divided into one of track 
control and of range indication, using a great-circular path between the 
points of departure and destination. 

The response of the track control system was examined with mecha- 
nization equations of as high an order as the fifth, and was found to improve 
as the order of the mechanization equation increased. Response to impulse 
type wind acceleration disturbances was satisfactory, but long-period sinu- 
soidally varying winds caused excessive errors. It appears probable that 
additional feedback loops, in conjunction with a mechanization equation of 
the fifth or higher order, will solve the track control problem satisfac- 
torily. 

In the range indication system, it was found to be impossible to re- 
move any forcing function terms of higher order than the time rate of 
change of acceleration. It appears probable that the range indication sys- 
tem, as postulated, will prove to be satisfactory, especially if the sensi- 
tivities are made variable, as a function of the angle of the input pendulum. 
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CHAPTER I 



INTRODUCTION 



The importance of developing a long-range self-contained system of 
automatic, mid -course, navigation for the control of guided missiles can- 
not be overestimated (Fig 1-1). As long range flights of piloted aircraft 
over enemy territory become increasingly hazardous, because of the neces- 
sity of travelling without fighter protection toward heavily defended cities 
provided with elaborate warning nets, automatic navigation can provide a 
method of delivering bombs without risking the lives of highly trained 
pilots. 

With the development of long-range jet-propelled missiles, automatic 
navigation should eventually permit the United States to bomb any place in 
the world from bases located within its territorial limits. 

For such a system of automatic navigation to be successful, it must 
have a degree of accuracy comparable to that achieved by human pilots. It 
must not be easily susceptible to "jamming" by enemy action. It must op- 
erate successfully at supersonic velocities, and in the upper reaches of the 
earth's atmosphere. 

At present, the problems of automatic long-range navigation are under 
consideration by several groups. Among these are: 

1. Baird Associates, Inc. 

2. Hughes Aircraft Co. 

3. Instrumentation Laboratory, M.I.T. 

4. Kollsman Instrument Company 

5. North American Aviation, Inc. 

6. Northrup Aircraft, Inc. 
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FIG. I- 1 RELATIONSHIP BETWEEN PHASES OP A LONG RANGE MISSILE T R A-iE.CTORV . 
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7. Rand Corporation 

8. U. S. Army Research and Development 
Sub-office (rocket), Fort Bliss, Texas 

Little work, however, appears to have been done on the derivation of 
theoretically optimum control equations or operating parameters. For 
this reason the writers determined to examine a long-range automatic 
navigation system now under development at the Instrumentation Labora- 
tory, M.I.T., as a means for exploring theoretically some of the general 
problems of automatic navigation. The Instrumentation Laboratory proj- 
ect, which is being executed under USAF Contract W33 -03 8ac- 13 969, is 
designated "An Automatic Navigation System - Project Febe," and will be 
the subject of a report soon to be published. This project has for its pri- 
mary purposes to determine: 

1. the feasibility of long-range automatic navigational 
guidance of bomber airplanes 

2. useful design parameters for a serviceable system 
for military use. 

Considerations other than military and tactical indicated that the sys- 
tem should employ solar tracking, a magnetic azimuth system, and a con- 
stant ground speed. The system is allowed to make no contact with the 
earth. Altitude is determined by the use of a barometric altimeter; the 
only inputs to the system are celestial observation, the magnetic field of 
the earth, and the various accelerations experienced by the airplane. 

This system, which has been installed in an Air Force B-29 bomber, 
causes the airplane to fly a definitely programmed great circle course at 
constant ground speed. Vertical accelerations are assumed to be of such 
small importance that they can be neglected. Figure 1-2 shows a functional 
diagram of the Febe automatic navigational system. During the period from 
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FIG. 1-2 

FUNCTIONAL DIAGRAM OF FEBE 
SYSTEM FOR AUTOMATIC NAVIGATION 
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June to August of 1948, the writers of this thesis were concerned with the 
Febe project in an under -instruction status. Through the courtesy of Dr. 
C. S. Draper, they were also permitted to attend a Seminar on Automatic 
(Celestial and Inertial) Long-Range Guidance Systems conducted at the 
Massachusetts Institute of Technology by the Scientific Advisory Board 
from February 1 to February 3, 1949. 

This background naturally led the writers to concern themselves with 
the form which a practicable service guidance system might take. The 
Febe system presents several technical and military problems, which 
arise principally in connection with the following: 

1. Constant ground speed 

2. The necessity of solar tracking 

3. Magnetic azimuth input 

4. Weight and size of components. 

Since a successful automatic navigational system is ultimately des- 
tined to fly in a long-range high-speed missile, it must be assumed that 
the weapon will be operating at close to the maximum range permitted by 
its size and fuel capacity. The navigational system, therefore, should not 
seriously lower the fuel economy. This suggests a system which will fly 
at constant, or nearly constant, airspeed. Furthermore, if the power plant 
consists of an athodyd, the airspeed must be maintained very nearly con- 
stant by the fundamental limitations of this type of propulsion. 

The limitation of the Febe system to a constant ground speed was 
largely dictated by the requirements of the azimuth system. The magnetic 
input could be eliminated through the use of two star -tracking telescopes. 
The problems of celestial tracking become increasingly severe, however, 
as missile speed is increased. At supersonic and near-sonic speeds, it is 
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no longer feasible to have an astrodome, in which to house the tracking 
unit, projecting from the fuselage. Thus, tracking through a flat window, 
with the attendant difficulties which this entails, becomes a necessity. 

Also, at high speeds, thick boundary layers and intense heating will exist 
along the surface of the missile, greatly complicating optical problems. 
Finally, the inclusion of a celestial tracking system materially increases 
the space and weight required by the navigational system — and weight and 
space factors become increasingly important with increasing range and 
bomb load of an aircraft. It appears, however, that approximately a thirty- 
fold decrease in the uncertainty levels of existing gyros would permit the 
maintenance of an inertial coordinate reference system within the missile 
through the use of such gyros alone, thus eliminating the use of celestial 
tracking. 

Great progress is being made in the improvement of gyros. The work, 
for example, that is being done by the group at the Massachusetts Institute 
of Technology under Dr. C. S. Draper shows promise of obtaining the nec- 
essary accuracy within the next few years. 

The model of the Febe automatic navigational system presently in- 
stalled in a B-29 is very bulky, and weights 1917 lbs. (See Fig 1-3). 

It practically fills the after pressurized compartment of the airplane. Of 
course, it must be remembered that no effort was made to decrease the 
size of this equipment, or to keep its weight at a minimum. Nevertheless, 
a missile held to a constant airspeed, rather than to a programmed ground 
spped, could, by using gimbal solvers, greatly simplify the elaborate trigo- 
nometric and speed computers which Febe requires. Much of the work of 
the computation systems might, of course, be performed on the ground 
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prior to the start of the flight, and be fed to the system through tapes, but 
this would introduce further complications. 

From all of these considerations, a self-contained inertial -gravita- 
tional system flying a great circle course at nearly constant airspeed 
appears to provide one of the most satisfactory solutions, if such a system 
can be implemented. 

After this examination of the Febe system, there remain, then, for 
use in the purely inertial system, a gyro -stabilized inertial platform, a 
clock mechanism to remove the earth's diurnal rotation, and two mutually 
perpendicular single -degree -of-freedom pendulous accelerometer units. 

In addition, there are gimbals and servomechanisms isolating the inertial 
platform from the motions of the missile containing it, and orienting the 
platform so that its axis is parallel to the polar axis of the earth, so that 
a reference vertical parallel to that of the point of departure (or any other 
convenient reference direction) is maintained independent of the motion of 
the missile. (See Figs 1-4, 1-5 and 1-6.) 

The three integrating gyros which are used to maintain a gimbal- 
mounted inertial platform "fixed" in inertial space, parallel to the earth's 
equatorial plane, are mounted so that they detect motion of the platform 
with respect to inertial space about the polar axis of the earth, and about 
two mutually perpendicular axes parallel to the plane of the equator (Figs. 
1-7, I~8). Servomechanisms, controlled by these gyros, rotate the gimbals 
so that the inertial platform remains fixed in inertial space. The earth 
platform, mounted on the inertial platform with the Identical polar axis, is 
rotated about that axis once in twenty-four hours, to transform the inertial 
reference system of the inertial platform into an earth reference system. 

3 
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FIG. 1-4. 

PHOTOGRAPH OF INERTIAL REFERENCE AUTOMATIC 
NAVIGATION SYSTEM MODEL. 
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FIG. 1-7. 

PHOTOGRAPHIC VIEW OF INERTIAL REFERENCE AUTOMATIC 
NAVIGATION SYSTEM MODEL WITH POLAR AXIS TIPPED AWAY 
FROM THE NORMAL DIRECTION TO THE BASE. 
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PHOTOGRAPHIC TOP VIEW OF INERTIAL REFERENCE AUTOMATIC 
NAVIGATION SYSTEM MODEL WITH CONTROLLED MEMBER GIMBAL 

SYSTEM REMOVED. 




Upon the earth platform, in turn, the controlled member is mounted 
in gimbals (Fig 1-9). The axis of the outer gimbal of this member is 
made to lie in the plane of the great circle connecting the points of de- 
parture and destination on the surface of the earth, considered as a 
sphere.* Then, about the inner controlled -member gimbal axis, the con- 
trolled -member is moved until it lies parallel to the surface of the earth 
at the point of departure. The inner controlled-member axis is placed 
perpendicular to the programmed great circle path so that, as the missile 
moves along the programmed path, the controlled member, by rotation 
about the inner axis alone, can be maintained at all times parallel to the 
surface of the earth below the missile. The angle through which the con- 
trolled member has rotated will then be a direct reading of the distance 
travelled by the missile along the programmed great circle. At some pre- 
determined angle between the controlled member and the gimbal connect- 
ing it to the earth platform, the missile may be said to have arrived at its 
destination. 

The two pendulous units are mounted upon the controlled member, one 
with its plane of motion in that of the programmed great circle, the other 
with its plane of motion perpendicular to that great circle plane. The lon- 
gitudinal pendulum, i.e. the unit with its input plane parallel to the plane 
of the programmed track, sensing longitudinal accelerations, produces an 
output signal which is then modified through the implementation of the lon- 
gitudinal mechanisation equation. The resultant signal is used as the input 

* This orientation is accomplished by rotating the controlled member on 
the earth platform to the correct meridian angle of the great circle track 
at the point of departure, L A Z (G.C.)(dep)^ and by mtin Z the outer gimbal 

axis to the correct polar angle (PA) / .. . of the great circle track plane 

with respect to the polar axis of the ^ earth. 
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FIG. 1-9 

PHOTOGRAPH OF CONTROLLED MEMBER MODEL WITH PENDULOUS UNITS IN PLACE 
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to a servodrive system that rotates the controlled member so that it 
tracks the local vertical. Similarly, the track control pendulum, sensing 
lateral accelerations; (i.e. accelerations perpendicular to the great circle 
track plane), produces an output which is modified through the implemen- 
tation of the lateral mechanization equation, and the resultant is used to 
position the missile control surfaces so that the missile remains close to 
the programmed great circle track. The altitude of the missile above the 
surface of the earth is maintained by a pressure type altimeter that sup- 
plies the essential input to the altitude control system. 

A Coriolis acceleration computer is used to compensate for the ef- 
fects of the Coriolis acceleration upon the action of the system. Geodesic 
acceleration is in general small enough to be ignored. It will be the pur- 
pose of this thesis to discover whether it is theoretically possible to main- 
tain the errors between the indicated and actual local verticals small 
enough to permit sufficiently accurate indication of the local vertical for 
the long-range bombing and missile guidance problems. Various mechan- 
ization equations will be examined, to see which offer the greatest chances 
of success. The most promising of these will be furthur examined to de- 
duce the optimum control parameters. 
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CHAPTER H 
KINEMATICS 



1 . Purpose 

It is the purpose of this chapter to derive the fundamental kinematic 
equations describing the motions of the missile postulated in Chapter I. The 
principle geometrical relationships are illustrated in Fig H-l , and are shown 
symbolically in Fig n~2. 

The following simplifying assumptions are made in deriving the kine- 
matic equations: 

a. Lateral and longitudinal motions decoupled 

b. Constant altitude flight path 

c. Spherical and homogeneous earth 

d. Controlled member located at center of gravity 
of missile 

e. Zero angle of attack and side slip angle 

The decoupling of the lateral and longitudinal motions divides the prob- 
lem into one of control of the missile in the direction perpendicular to the 
programmed great circle and one of indication of the distance travelled in 
the programmed great circle. The kinematic equations are therefore de- 
rived separately for lateral and longitudinal motions. 

Since the flight path is to be considered as having constant altitude above 
the surface of the earth, the only vertical component of acceleration is the in- 
ertia reaction acceleration of gravity. With the earth assumed to be spherical 
and homogeneous, this componentis a constant over the surface of the earth, 
with the numerical value of 115,920 feet per minute per minute (equivalent to 
32.2 ft/sec ). The assumption of a spherical earth also eliminates any geo- 
desic acceleration terms from the kinematic equations (ref par 1 App A). 
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FI6.II-1. GEOMETRICAL RELATIONSHIPS AMONG PHYSICAL QUANTIT 
ASSOCIATED WITH TRACKING IN AN INERTIAL GUI OANCE SYSTEM. 
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A Coriolis acceleration computer is included in the fundamental guidance 
system (Fig 1-5), and is assumed to work ideally (ref par 2 App A ). The 
kinematic equations, therefore, do not contain the Coriolis acceleration. 

With the controlled member located at the center of gravity of the mis- 
sile, the gimbal isolation system can be considered to remove entirely all 
effects of the roll and pitch of the missile. No accelerations reach the con- 
trolled member as a result of angular accelerations of the missile about its 
body axes. 

Considering both the angle of attack and the angle of side slip as zero, 
the velocity vector of the missile always lies along the longitudinal body 
axis of the missile; that is, V[( a i r ) - 14 ] lies along X M (see Fig n-3). This 
is taken as the control direction (CD). The kinematic equations for lateral 
motion are derived first. 

2. Kinematic Equations of Motion Perpendicular to Great Circle Track * 

The geometrical elements of the simplified lateral guidance problem are 
illustrated in Fig II -4. Figure II -5 defines the directions and angles that are 
used in the derivation of the kinematic equations. 

As can be seen from Fig II-5, the following relations are true: 

(C)A( Vert x tC ) - (DC)^ true >(tc) “ (DC)^ n d)(tc) ^ ^ 

• • 

= ^E (C^A(v er f)(tc) (n-2) 

tan(DC)( trae)(tc > ) (n-3) 

♦Irhis development follows that of John Hutztnlaub, Notebook, July, 1947. 

20 







n 



FIG.n-3 MISSILE COORDINATES 
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FIG.n-4 GEOMETRICAL ELEMENTS OF TRACK CONTROL PROBLEM. 
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i • 

R E (C)A (Vcrt)(tc) = - g iR tan (DO (true)Ctc) 
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But, from (II- 1) 

tan (DC^ (t rL1 e Ktc) = tan [(C)A^ Verl){tc) + (^C) (ind)(tc) ] 

tan(C)A(Vcrt)(tc> + tan (DC) ( ^ nc |^ c ^ 

1- tan (C) A( Ve rt)(to tan(DC) an d K tc> 

(II-6) 



ton (OA^ Vert ^ tc) 

l-tan(C)A (Ver t)(tc) tan (DC) (cnd)( tC ) 

tan ( DC) (indutc) 

1- tan (C ) A (v/ert) (tc) tan (DC) ((ndKtc) 

(U-7) 

This is a theoretically correct expression, if no vertical accelerations are 
present, or if vertical accelerations are assumed to cause a variation in 
g^. However, if the size of (C)A^ vert ^ tc j does not become greater than 10 
milliradians, an error of no more than one part in 30,000 is introduced when 
(C)A^ ver ^ is substituted for tan (C)A^ yer ^^ c ). This should be true for all 
cases of practical interest. If this is done, eq (n-7) becomes: 



Therefore, 

R 



lC)A(vertktc) 



IR 



CO NilUi- f: ' i 

J i^U 



lC)A (Wcrt)( tc)+(-^- ) 



(C) A 



CVertXtc) 



R e l-(C)A( Vert ) t t c) tan ( DC) (ind)(tc> 



v r f > l- 



tan (DC) (i . n dutc) 



(C)A (vertH^cj tan (DC)^ n< j|j(t. C ) 



grp x 

II p— - is defined as W^, then eq (II”8) becomes: 

/r\A + W a (C)A (V ertHtc> 

(CMWrt nE l-(C)A (Vert ,(tc ton(DC) (lndKtc , 



- w 



tan(DC) (lnd|(tc) 



nE l'(C)A (VertH tc> ton (DC^( Lnd ^(^ C ) 



(n-8) 



(« -9) 



If (DC)^ in( j^ te ) is smaller than fifteen degrees,* ^-')(i n d)(tc) can ^ 
substituted for tan (DC) * or order theoretical solutions. Also, 
(C)A^ yert ^ tc ^ tan (DC ) * as die product of two small angles, becomes 
very much less than unity, and can be neglected. Then: 



(C)A 



(Vert)ttc) 



+ < E ©A 



(VertHtc) 



= W nE (DC) (ind)(tc) 



(n-io) 



Using the same small angle assumptions, and the additional assumption 
that Ar 2 ( cm j remains less than fifteen degrees, eq (11*4) can be written in 
the approximate form: 



lE-(cm)] (tc) “ ^[E-(air>] (tc) ^ 



[(air)-M] A Z.(cm) 



(n-u) 



* Ref par 3 App A 
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Equations (n _ 2), (n-10) and (n ”11) are the basic kinematic equations that 
are used in all the derivations of Chapter in. 



3, Longitudinal Motion Kinematic Equations * 

The longitudinal problem is one of range indication. The geometrical 
elements of the simplified range indication problem are illustrated in Fig 
II-6. Figure n-7 defines the directions and angles that are used in the der- 
ivation of the kinematic equations. 

As can be seen from Fig n-7, the following relations are true: 



(DC) (t ru .e)(Jortcp = (C)A(vert)Uong) + lDC)(ind)(lon^) 



LE- (cm)l(lon^) 



k r(true) 



(n-12) 



(n-13) 



tan(DC)( true ^| on( ^yr 



a 



(E-(cm)] ( loncj) 

9ir 



(n-14) 



Taking the tangent of eq (n - 12) 

tan <DC) (true >< lonq) = _^ Qn(C)A Cvert)(lon« } )- t -ton(DC) (CndHlon<a , 

(trucKlo.*) 1-tan tC) A (vert)Uonij) tftn(DC) (llK j M | 0 ^ ) 



(H-15) 



Equation (II“15) is a theoretically correct expression but, as in the 
lateral case, all practical conditions find (C)A^ ver ^j 0 ^ smaller than ten 
miUiradians. Then, eq (n-15) simplifies essentially to: 



♦This development follows that of John Hutzenlaub, Notebook, July, 1947. 
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FIG.IT-G GEOMETRICAL ELEMENTS OF RANGE INDICATION PROBLEM. 




} 




$ 

4 

fig. n-i 

LONGITUDINAL INDICATION DIAGRAM 
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tan (DC) 



^^^CVert)Clong'> 
(truOdon^, - l- (C ) A(Vert)(lons) 



* tan (DO^^Ioh^ 
tan (DC)^ n< ^ )( | on< ^) 



(n- 16 ) 



From eqs (n-13) and (n-14), 

• • 

-•» / rv rt\ _ A r (true)^E- _ 

tan(DC) ctrueH , on< ^ - g 



IR 



(true.) 



Substituting eq (II-17) into eq (II"16) 

• • 

^r(true> ^ ^ (Vert) ( long) * tan long) 

W n £ ^"^^(Vcrt) ( loncj) tan ( DC ) (Lnd)(| 0 n<j) 



(n-i7) 



(11-18) 



From this there results, as the general longitudinal kinematic equation, 
V _ (C)A (VertU | onq) 

r(true) nE ]-(C)A( Vcrt ^ (( on ^) tan (DC) ( £ nc ^ (| 0 n<j) 




ton (P£) (tnd)( long) 

1 - (C) A( Vcrt ^ ( 

lont} 1 ) tontDC) (indu 

lonq) 



(n-i9) 



When the same small angle assumptions are made for the longitudinal 
kinematic relationship that were made in writing eq (n-10) in the lateral case, 
eq (n-19) reduces to 



ArCtrue) - W„ E (C)A( V€rtU | on ^)+ W* £ (DC)(ind>(lon0^ 



(n- 20 ) 
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This equation and the relationships 

hr(true) = A r (ind) ~ (^ A( Ver t)( |on<j) 



(n-21) 



and 



A r(a pp) = A r(Ln d)^ (DC) (tndHlon^) (11-22) 



form the basic kinematic equations that are used in all the derivations of 
Chapter IV. 



30 




CHAPTER III 



THE TRACK CONTROL PROBLEM 



1. Introduction 

The simplified kinematic equations developed in Chapter II 
for the track control problem are restated here for convenience. 



Using these simplified kinematic equations, various mech- 
anisation equations are investigated in an effort to determine a 
physically realisable and reasonably accurate mechanisation of 
the track control problem. As previously stated, for the purposes 
of this study, the track control problem is considered as de- 
coupled from the longitudinal problemil)^ and vertical accelera- 
tions are neglected(2) . it is further assumed that the missile 
operates with perfect control surface servomechanisms and that 

1 - ref par 2 App A 

2 - ref par 4 Aj>p A 




(m-i) 



(C)Arvertxtc) +w u \ (C)A 



(ve»-t)(tc) 





(m-3) 
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0 




iv'P' 




l 

L. 



(3) 

it possesses perfect aerodynamic response' ’ . Finally, the pendulous 

accelerometers are considered to indicate instantaneously the direction 

(4) 

of the resultant acceleration' ’ . 

The purpose of the track control mechanization is to reduce toward 
zero any angle (C) [(Vert)^^ ](t c ) that has been introduced by lateral 
motion of the missile away from the prescribed great circle track, and 
to maintain this angle sufficiently small to insure arrival within one 
mile of the destination. The control is achieved through operations upon 
signals received from the measurement of the available angle 
(C) [(Vert)app]^^, the pendulum angle, in order to correct the head- 
tag angle, A [ X(cm) - (CD)]' 

Throughout the derivations which follow, it is convenient to refer 
to Fig n-5, which indicates the relationships among the important angles. 



2. Proportional Control 

The simplest possible mechanization is proportional control, in 
which the heading of the missile is changed in linear relation to the 
pendulum angle as measured between the pendulum and the controlled 
member. This case is therefore considered first. The mechanization 
equation for this concept can be expressed as: 




ZCcrrD 




LC(CD)][AAJ 






(m-4) 



(3) ref par 5 App A 

(4) ref par 6 App A 
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Using this mechanization equation, and the simplified kinematic 
equations, numbered one through three, a performance equation 
for this system is now derived. 



From eq (III-3), 



Vre.^raftc) ^CE-fcm)J (tc) 



iz(crn) 






M) 



(m-5) 



Integrating eq (III-l), assuming the constant of integration as 
zero^D, and substituting into eq (III-5)> 



A =. £ - (a.\ y)3 (t c) 

l(cm) 



M C)A (Ver 



t) (tc) 



(ni-6) 



Cl' dir)- MJ 



From eq (III-2), 

% 



fD('\md)(Vc) 



_ (c)a 



fV'rtKtc) 



w 1 



"HO A 



(Vcrt)(tc) 



(ni-7) 



1 - par 7 App A 
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Substituting (Ill-b) and (III-7) into the mechanization eq (II 1—4 ) , 



V 



C£-fair)](-fcc) *-£ 



-z (C)A 



V 



CrAir)-M3 



(DAj 

W‘ 

V V Nl 



CvertUta — c J ^ // ~"Vvert)(tc) I A 

^CcfcoUCAA) 1 \ /i T \L;M( Vert )^c) 



(lll-b) 



From this, the performance equation can be written: 



(C)A 



(Vert) fire) 



e E w: 



V 5 

V f(dir)-M] ^[C(CD)]tAAl 



(C)A, Vert)ftc) + W; e (C)A (vert)(tc) 



w. 



i 

fH E 



V 



Le - (air)] (id 



5 



CcrCD]J[AA] V ! 7 air}-M) 



V r 



(I ll- 9 ) 



Although this equation provides a stable control system, a 
forced error results from the lateral component of the velocity 
of the wind, V ^ _ (air)] ^ y This, obviously, is unsatisfactory 

for the control of guided missiles, because no large steady-state 
errors can be tolerated. 



3, Integral Control 

Next, let us examine a system that changes the heading of 
the missile in accordance with the integral of the pendulum angle. 
The equation for this mechanization can be written 




(II 1-10) 







‘^U(w t .ML 



Differentiating (III— 6) , and substituting with (III-7) into 
(III-10) differentiated, 



V[E-fair)]ftc)“ -RelOAfc 



(Vert) (tc) 



V ' 



= 5 



[cfc D^AA] 



[(air)-M] 



^*^A(Vert)(tc) 

~wT 



+ (C)A^ r 



tXtoj 



(111-11) 



or 



(0A (Vert)(tc) -h 



w l 

Y y ne 



fc)A 



W z 

WE 



(VerO(tc) 



V r 



IR- 



S . V 

v - / CcfcD)3CAA] v C^a.ir)-M] 



- ~h 1 



g _I_C .V 'Ce-^rOftc) ( | | | -1 2 ) 

■° ' ^CcCcd)3Caa] Y C(Ojr)-M] 



I It 



This system, although the forced error is now due to the 
acceleration of the wind, has no damping term. S[(q)(cd)] [AA] 
can be made large to reduce the size of the forced error, but 
this results in an undamped period of oscillation of approxi- 
mately 84 minutes for the controlled member. Any residual forced 
error (since it is physically impossible to make S[(c)(CD)] [AA] 

infinite), causes this oscillation to appear. 



4. Second Integral Control 

In order to reduce the forced error, consider next a system 
that changes the missile heading in obedience to the second inte 
gral of the pendulum angle. Here, the mechanization equation is: 




dt dt 






Cl l 1 - 13 ) 



ae 



Differential ing (III —6 ) twice, and substituting with (III-7) 
into (III-13), differentiated twice, 



-(V>r)y tc) 



— R. (C)A 



(Vert)ftc) 



fC)A rvert )(tc) 



V 



[to.‘i r)- m] 



tcfco)^ C^a3 



W MI 



+ fc)Acv.rt)a 



cj 



(111-14) 



or, 



(C)A 



+5 



V, 



(ver-t) Ctc) > U[c(cd)3Caa] 



CCdJr)-M] 



te 



(Ok 



(vert)(tc) 



f 5 



V, 



U*,r)-hG /(-U W 

CCCcdD-C aa} r-> V'- / /^A(vertKt^ v < 



R 



R Ce ' ft,r)5ftc ) (in- 15 ) 



By Routh f s stability criteriafD, this is unconditionally 
unstable . 

5, Derivative Control 

In order to give a complete presentation of possible types 
of mechanization, consider the heading angle to be controlled in 
accordance with the derivative of the pendulum angle: 




( I I I -16) 



1 - Routh, " Advanced Rigid Dynamics'’. 



Dif i ere ntia ting (III-7), and substituting, with (IXI-6) into 
(III-16), differentiated, 



V, 



[e-cairn etc) 



-r.toa 



(Vert)ftc) 



V 



= 5 



Cfeir)-K] 



[cCco)3CAA] 



roA 



w.) 



(Vert)Atc) 



h-(c)a 



(Vert)(-tc) 



NE 



(111-17) 



From ( III— 17 ) 



^)A lV e^c) + W, 




R, 



v 

A V 1 






V 



[e - fair-)] (X c ) 



CcfcD)][AA) v Cte>r)-MLl < 



V, 



[fair)- Mj 



T 



5 

v 'rv-:'co)][AAi 



( 1 1 1 -Id) 



By Routh’s stability criteria for a cubic, this is uncondi- 
tionally unstable. Moreover, there is a forced error due to the 
velocity of the wind. 



6. Pr oportional Plus Integral Control 

As shown in paragraph 2, proportional control produces a 
solution which, although it has a forced error caused by wind 
velocity, is nevertheless stable, while integral control (para- 
graph 3 ) develops a solution that does not have any forced error 
caused by wind, but which also has no damping, and is therefore on 
the verge of instability. Thus it might be expected that by com- 
bining both equations, a stable system without forced error due to 
velocity of wind might result* The resulting mechanization equation is: 




(Ml-19) 







Differentiating eqs ( III— 6 ) and (III-7) and substituting with 
eqs (III-7) into (III-19) differentiated. 



• •• 

VcE-kirl] (tc) ~ Re (C)A (Vert^tc) 



V 



= 5 



ffC)A 






[CfcD)} fA Al 



fvertV-fcc) 



W 1 

V * N£ 



+ (C) A 



<V«rt)(tc) 



+ 5 



C£ (c d)J Taa3 



(C) A( Vfrt)ftc) .s 



(111-20) 



From eq (III-20), 






w z 

VV NE 



R 



Cc (cd")J [aaJ 



V, 



"h 



Ccfcofl La a}! 



C^i r)-M] 



w: 



rc)A 



('Vert ) (-fcc) 



+W Me (C)A 



/Vert)(tc) 



~hW we 



^CC(CD)] CAAl Yc)a 

5 £c(cD)] [AA] 



(Vert)rtc)' 



w 2 

VV Nit 



CC(CD)] [aaD V c ^»V)-m] 



V 



C (111-21) 



By Routh* s criteria of stability, this system is stable if: 
to SfftVfcoa caa] aM<1 5 twfcearA ij *r«. positive ^u^bers. 



Cb) 



w l 

vv ne 



s v 

*1*0(0 C A A] R’a.rrVM'} 



■>o. 



Furthermore, as expected, this system has no forced error 
resulting from wind velocity, although it has one caused by the 
acceleration of the wind. This suggests adding second - integral 
to the system just discussed. 



* 



88 



<-0 Worn :al 



7. Proportional Plus First and Second Integral Control 



A] ^OciwlKtc) S R(co)] [A Aj / ^ ^ fc C* AJ 



Substituting eqs (III-6) and (III-7) differentiated twice, 
eq (III-7) differentiated, and eq (III-7) into eq (III-22) 
differentiated twice, 




M 

V 



[g- (Airl] (-kc) 



- P. (C ) A 



(Vcrt)(tc) 



V 



=5 , ••• 

^[c^co^DLA / 



A3 



CC&'iO-mJ 



y MM 

' ' 'Yvcrt) ftc) | 

ME 



(\ r\d)(t<0 



dt dt 



(Vert) (*c) 



^5 



Cc (cd)j [a a 3 



(C)A 



(Vert) c) 



w z 

we 



+^)A rw ^ fc) r-H5 



(Ver't) ftc)| '-'£(co) 3 CA A3 



fc)A 



(Vert) (tc) 



W 1 

ME 



+■ (OA (vc 



#-t)(tc) 






( I 11-22) 



(III -23) 



From eq (III-23) 



(C)/\ Ve ,^tc)+ 



IK- 



S ..V 

^CC(C0)3lA AO V C<^« r )“M3 



4- 



^ Cc (cd)3Ca A3 



5 



[c (c d) 3 £& A 3 



(C)A (vtrt)fe0 -H W* G + 



5 



Cc(co)3CAA3 



5 



(c(co)J Ca A 3 J 



(0 A 



(v«,t)ftcr 'c 



Cc(CP)J CA A 3 



w;« (c)A 



(Vert)(tc) 



£c(cp) 3 Ca a J 



s W 1 

I ^CcCcp^ Caa 3 y / l (C\A — ^ 

‘ c W Ne , 

Orw. -n rV T t x rv dr.'ll Cl A *1 V r 



V [E -r 4 , r)J (1 1 1-2+) 



Cc(cD)J CA Aj 



Cc (cd )3 ca a 3 v cr Al r )- m 3 



According to Routh’s stability criteria for a quartic, this 
system is stable if: 

(a) The three sensitivities have the same sign. 
(They may be all positive or all negative), 



(b) 



W 5 • 

VV ME ^CcfcC» 3 CA AJ ^ 



)K. 



'fcfeoj CaaJ 
£c(cd) La aj I 



"b 



(cc >)3 CA Aj 

y 



>o 
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In this system no effect of the wind of lower order than the 
second derivative of the wind velocity enters into the forced 
error. It seems unlikely that a forced error of this nature could 
have a magnitude large enough to be of practical importance. How- 
ever, if this error should prove important, it is possible to go 
to higher order performance equations. One further case is con- 
sidered. 



8, Proportional Plus First, Second and Third Integral Control 




dt dt— 5 




dt dt dt 



Substituting eqs (III-6)and (III-7) differentiated three 
times, eq (III-7) differentiated twice, eq (III-7) differentiated, 
and eq (III-7) into eq ( III— 25 ) differentiated three times, 



(MI-25) 






V 



Ce -foV)3[tc] 



-R„(C)A 



(VertKtc) __ 



V. 






-5 






CcfcpntA Al 




+rc)A 







(C)/\v.h 



CcfcpflCA A # J 



w 1 

M£ 



<v --H r C)A 



, -1-5 

A/ert)^0| ^Cc(co)J [A A] 






94 



w 



l-i-Q 

i ' VWM <Wytc)( ' 



'CcfcD/J [a’aJ 






we 



W 

T T N£ 






(III -26) 



From eq ( III— 26 ) 



(0 



IK. 









+ V’* 1 ™ | c)Aw<w+lw:+ Sssasia 



s> 



:cCcp)3C‘a AO 



^CefcoOO« 






w 1 

’▼Ml 



s 

5 



CcCcd)} Ca ‘aJ 



-h 



s 



Cc.(co)3 Ca A3 



) 



Wco)J t*A A 3 



Tcrco)3L'A‘ *X3 



^A^ rtV+c) -HW Ne 



* s 



{cfcoDTA aJ 



5t 



c(coY}Ta a/ 



(w;. f- ° ,KAM 



► /If) A ^ Wwe 

f lW/A rVcrt y t0 _ r ^ 



TcfcoOtA A]/ 






V, 



(in-27: 



S ' \/ v [e-^irO 

C^cdDCa aJ V>V)-mJ 
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By Routh f s stability criteria for a quintic, this perform- 
ance equation is stable if: 

(a) All coefficients of the homogeneous 
equation are positive, 



(b) 



5 

5 , 



tCtcPlKAA] 



w 



X 

ME 



cc(coBr a ad 



s -l 

^CciCD) 3 rAAl ' 




w 



X 

Ne 



Cc(cd) 3 Ca O 



<3 



IR 



S V 

'■''frfc. pUCaa'3 v [(WVm3 



+ 



s 



(c(co)3 Ca a 3 



S 
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1 R. 



+ 



s 



Co(cD)KAA 3 



^Cc(cd) 3 C A A3 r)-M3 CA A3 




Cc(cp)JCa a3 



Cc(cd)] C a a3< 



r 



9 



I R. 



J Ac d33 C A a^X ( iir)-M3 






Cc(cd)D CA a] 



5 



Cc(cd)3Ca a], 
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1 

Ni 6 
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• • ♦ • • . 
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S 



rc(CD)]0 A 3 



[c(co>3C'a* a*3 



Cc(cd)3 Ca a 3 





C c (cd)JC a 



Ajl 



s s 7" 

x/ 1 nJ Cc(cp)3Ca a 3 | ^cc(cp)][a a3 

> ^0Kco)]CaaJ A-’-coucaaV 



Me c 

^[cfcotfCA a’ 3 / 



IR- 



H” 



S • x/ 

^Tc(CP33LA A] 



Cc(cd>3La A 



V S 



[cfco)3CA A 3 



ccfconcAAi 

v " e 5 

Cc(co)]C'a A3 




;c(com'A A3 



5 



Cc(cd)3Ca a3 



NE 5 5 

lc(cd) 3C A A 3 ^CContA A3 V 



These restrictions permit the sensitivities to be either 
positive or negative, provided that they all have the same sign. 

It is of interest to notice that in each stable performance 
equation, the forced error has the same type of coefficient, with 
the order of the wind derivative increasing from velocity for 
proportional control to any derivative desired, as the perform- 
ance equation increases in complexity. This coefficient is of 
the form: 



[C(C 0 )HA aJ 



Cc<cd> 3 Ca a 3 
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The forced error is thus seen to decrease as the missile airspeed 
increases, and as the sensitivity S[C{CD)] [M] incre ases. Theo- 
retically, if this sensitivity increased to infinity, there would 
no error, either transient or in the steady state, caused by 
the action of the wind. This is similar to the conditions found 
with integral control, eq (III-12), where, as the sensitivity 

S [C ( CD ) ] [AA] increases to infinity, both steady state and trans- 
ient wind errors reduce to zero. 



— Track Control Corre ction as a Function of w<nH 

Since the mechanization eq (111-22) for proportional plus 
first and second integral control appears to provide a stable 
system without excessive forced error, this system is considered 
further. The performance equation for this system, eq (III-24), 
can be expressed in operational form: 
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If (C)A( vert j is expressed in minutes it is a measure, in 
nautical miles, of the transverse linear distance of the missile 
from the desired great circle track* The oresence of the opera- 
tional symbol in the numerator of the equation indicates that 
there is no error, in the steady state, from an accelerating 
wind* If it is desired to find the missile response, to 

an accelerating wind, it can be calculated if eq (111-2$) is 

A *7 f q rti ) 

multiplied by an expression for ' : The result 

TcTA (vert)(tc) • 



yields the heading of the missile as a function of wind velocity. 
Manipulating eqs (III-l), (III-2), (III-3), and (III-22), and 
expressing the result in operational form, 




Combining eqs (111-2$) and (III - 29 ) gives, 
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The negative sign shows that the heading Az( cm ) of the air- 

e * 

plane is correct to oppose the tendency of V [g _ (air)] to blow 
the missile off the great circle course. This is true even if 
the sensitivities are negative. 

10. Track Control Closed Loop System 

The relationship of eq (III-30) can be obtained from a study 
of the over-all lateral system as a servomechanism. In the direct 
servo closed-loop system, the angle (C)A( vert j does not appear, 
but it can easily be obtained, as indicated in the accompanying 
Fig III-l. 



In this figure: 



(PF)i has the value — ±— 

«IR 



(PF) 2 



^Z( cm) 

< DC >(ind)(tc) 



and is made up of the sum 
of two performance functions 



(PF) a 



SiMlitcj and 
( DC >(tru*)(tc) 



(PF) b 



^Z(cm) which is the mechanization 

(DO (ind) (tc) equation in operational form 



(PF)o, in the feedback path, is (- ^ [( a ^ r ) . z i H )p 
J SIR 

providing the proper function for the com- 
binor, which is physically the pendulous unit 
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FIG.nr-1. TRACK CONTROL 5YSTEM A 5 
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The pendulous unit, if (DC) ( trU e) remains smaller than about fif- 
teen degrees, solves the equation 
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(4rae)(tc) 



IR 
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( M 1 - 31 ) 



This equation is derived from eqs (III-l), (III-2), (III-3) 
differentiated, and the relationship that 

(DC) aru€ )( tc )"~ ^C)A (VertHtc.) (III -32) 



This last equation can readily be obtained from an inspection of 
Fig II-le 

Solving for the performance function (PF) a> using eqs (III-l), 
(III-2), ( III— 3 ) , and (111-32), 
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The mechanization equation, which in this example is eq (III-22) 
expressed in operational form, is 
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The product of eqs (III-33) and (III-34) gives 
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but, from the theory of servomechanisms, 
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Solving this equation gives eq (111-30), as before. The stability 

criteria for eq ( III— 30 ) are identical with those for eq (III-24). 

Since the sensitivities can have either sign, a necessary condi- 

glR 

^ [E - (air)] 

if the sensitivities are chosen as negative. If this is expressed 



tion for stability is that 
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1 - Brown and Campbell, "Principles of Servomechanisms 




then, for stability, 
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11, Numerical Solutions, Using Two Quadratics 

Using the relationships that have been established, and 
substituting numerical values for the coefficients, it is possi- 
ble to solve for the various angles as a function of the accelera- 



using a range of values for the sensitivities, it is possible to 
select the most satisfactory values, and to estimate the effec- 
tiveness of the system. This is now done. 

The homogeneous portion of the performance equation, (III-2/*) , 
is a biquadratic. It can then be written 



which factors the equation into two quadratics. Manipulation of 
the equation indicates that if the damping ratios and natural 
frequencies are varied, the higher damping term associates itself 
with the shorter of the two frequencies. This result then appears 
essentially as that of a lightly damped long period quadratic, 
which is undesirable. For this reason (DR)^ and (DR ) 2 » sad 
and W 2 have been chosen equal for all solutions of the equation, 



tion of the wind V j-g (air)l • From a sequence of such solutions, 




(IM-39) 
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and equations are now derived to facilitate the selection of 
values for the sensitivities to make these limitations hold. 
Expanding (III-39), in the general form, gives 

p 4 +[2 (DR)W+Z(DR\yjp 3 + [w:+W> Hip 1 + 

[z(D^y l w x +z(DR)w,w;] p+ w; w; =o (m . 

Wi W 2 

If - «-=- - (FR), and (DR) i - (DR) 2 - (DR) then, 

W NE W NE 

comparing the coefficients of eq ( III-40) with those of eq (III-24) 
gives the following relationships: 
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When these four equations are solved simultaneously for the 
sensitivities and for the damping ratio, the results are: 
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(III -45) 



(III -46) 



(111-47) 



(III -48) 



It is noticed that the choice of a damping ratio immediately 

fixes all sensitivities and the frequency ratio (or that the choice 

> 

of a frequency ratio, alternatively, fixes the damping ratio). As 
the frequency ratio increases; i.e., as the quadratics have pro- 
gressively higher natural frequencies, the damping ratio simul- 
taneously increases without limit. Conversely, as the frequency 
of the quadratics approaches that of the #4 minute period "earth 
pendulum", Wjjg , the damping reduces to zero. The sensitivities, 
at the same time, increase without limit. 
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12. Numerical Solution Using Four First Order Terms 

A second method of writing the quartic equation, which can 
be examined easily, and which may yield results of practical 
interest, is that in which the equation is reduced to four decay- 
ing exponentials. This equation is of the form 

(p + «y (p + (s) (p + «b) (p + <f>) = ° 



When expanded, this becomes 
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(ct) , (ct) z (cr) 3 (cr)^ w ( 111 - 50 ) 

The general relationships among the coefficients of eq ( III— 50 ) 
and eq (III-24) give 
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The special case which Is considered here is that for which 

% 

(CT)^ - (CT )2 * (CT)^ * (CT)/^ - (CT). Under these special circum- 
stances, the simultaneous solution of eqs (III-51, 52 , 53, and 54) 
yields 
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Evaluating (III-57), (111-58), and (III-59), using eq (111-56) 
with the smaller (CT), 

(111-60) 
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Cc(co>lC*X3 

This gives negative values to the sensitivities. The alterna- 
tive evaluation would make the sensitivities positive, but would 
result in characteristic times f$r longer than could be tolerated. 
The system, with negative sensitivities, is stable by Routh f s 
criteria. 

13 > Numerical Solution Using Quadratic and Two First Order Terms 

The final method of examining a quartic equation which may be 
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of interest is that in which the quartic ia broken up into one 
quadratic term, and two decaying exponentials, 

(p'+Z(DR)Wp+W‘)(p^)fp+ ^ oi.-fei) 



Upon expanding eq (III- 63 ) this becomes. 
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The general relationships among the coefficients of equations 
(III-64) and (III-24) are 
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These equations can be solved simultaneously for the sen- 
sitivities, and for the relation which connects (CT)]_ , (CT>2 
W and (DR). It is convenient to use the relationships: (FR) « 



W 



NE 



v [ (CT)PR 3 - 2w 

resulting equations are: 
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These are evaluated for the single case of (DR) - /£ , with 
[(CTJPR]! - [(CT)PR] 2 - 0,1. Then, from (III-69), (FR) . 0,307 . 

This gives : 



\ R. 



Kfco>uX ' M W NE 37.G4 



CC(cd)Ha'aJ 



9,. jgOS 

V, 



[air)-M] 



0.714 



'CcCc D)][A AJ 



V, 



[(Vir)-NV} 



(HI-73) 



(ill -74) 



(Ml-75) 



Notice that once again the sensitivities are negative. 

14 . Plots of Results 

It is now possible, using various coefficients for the 
equations of paragraph 11, or the indicated coefficients of the 
equations of paragraphs 12 and 13 , to obtain numerical answers 
by solving the equations with various wind inputs. This was done 
on the Rockefeller Analyzer at Massachusetts Institute of Technology. 

It is not necessary to solve for various missile speeds, since 
the speed of the missile can be combined with the heading angle of 
the missile A Z ( CM j. This is possible because of the linearisation 
assumption which was made in deriving the kinematic equations, and 
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hold* only if the missile airspeed is very much greater than the 
velocity of the wind. 

Plots were made of the angular error in minutes of arc of 
the position of the missile perpendicular to the great circle 
track. Considering the earth to be spherical, this is directly 
a measure of distance from the track in nautical miles. For 
convenience, on all plots, the unit of angle is the minute of 
arc, and of time, the minute of time. The foot has been taken 
as the unit of di stance , 

The first input of wind which was used in solving the equa- 
tions was a rectangular pulse of wind acceleration, or a constant 
acceleration of the velocity of the wind from sero to some finite 
value. The duration of the pulse was taken as one minute. Next, 
responses were taken for a step function input of wind accelera- 
tion. Finally, responses were determined for sinusoidal inputs of 
wind velocity using several frequencies. Plots of the results 
follow, using various sensitivities in eq (III-22): 
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NUMERICAL CONSTANTS USED IN OBTAINING 
RESPONSE OF TRACK CONTROL SYSTEM TO PULSE AND 
STEP OF WIND ACCELERATION AND CIRCULATING WIND 



TABLE m-1 
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FIGURE m-2 RESPONSE OF TRACK CONTROL SYSTEM TO A PULSE 

DISTURBANCE OF WIND ACCELERATION 
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FIGURE m-4 TRANSIENT RESPONSE OF TRACK CONTROL 
SYSTEM TO A STEP INPUT DISTURBANCE OF WIND ACCELERATION 
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WIND (20 MINUTE PERIOD) 

(FOURTH ORDER SYSTEM) 
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RESPONSE OF TRACK CONTROL SYSTEM TO A 31.4 mph CIRCULATING 

WIND (30 MINUTE PERIOD) 

(FOURTH ORDER SYSTEM) 
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CHAPTER IV 



THE RANGE INDICATION PROBLEM 



1, Introduction 

The simplified kinematic equations developed in Chapter II 
for the longitudinal problem are restated here for convenience: 



Using these simplified kinematic equations , various mech- 
anisation equations are investigated in this chapter in an effort 
to determine a physically realizable and reasonably accurate 
mechanisation of the range indication problem. The assumptions 
are identical with those for the lateral problem. The missile is 
assumed to move with instantaneous response to wind forcing 
functions. 

The purpose of the longitudinal mechanisation is to reduce 
toward aero any angle (C) [(Vert) (ind) (long)] which has been in- 
troduced by motion of the missile along the prescribed great 
circle track, and by measuring the angle between the reference 
vertical and the indicated vertical, to obtain a measure of the 




(IV-1) 




( tV-2) 




(IV-3) 
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location of the missile. The orientation of the controlled member 
is achieved by measuring the angle (DC ) ( i nd ) ( i on g ) between the 
position of the pendulum and the controlled member, and by oper- 
ating on this quantity to position the controlled member with 
respect to the reference vertical. The longitudinal reference is 
chosen as the vertical of the departure point, ( Vert ) (dep) • Since 
the value of (W) (ind) (long) is dependent upon the position of 
the controlled member, as well as of the pendulum, it is more 
convenient, in deriving mechanizations, to use the angle A r ( a pp) 
between the reference vertical and the indicated vertical of the 
controlled member, than to use (DC ) ( i nd ) ( long) directly. A r ( app ) 
is physically measurable, since it is the sum of (DC ) ( i nd ) ( long) 
and A r ( ind ) , which are both known. It has the advantage of being 
independent of the position of the control] ed member. 

Throughout the derivations that follow, it is convenient to 
refer to Figure II-6, which indicates the relationships among 
the important angles. 



2. Trial Mechanizations 

A. The first mechanization that is considered is that in 
which the second integral of A r ( a pp) moves the controlled member 
according to A r ( ind ) and its first two integrals. That any simpler 
case produces an unsatisfactory solution is indicated by the mech- 
anization equations for the track control problem. The mechaniza- 
tion equation is 
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[c(c.*r>V)[A A] 




^ f I ^ 



/ ■ * » I 

v <J I -j h 



> u ii' 






Differentiating twice, 



5 



[cfcm)][AA3 



A r find) • 5 



Lc(cm)] [aaJ 



A r ( m d) ~^5 



tc(c^T!CAA3 



• I 

A r (,- r ,a) = A f 



(&pp) 



(iv-5) 



Substituting ( IV— 3 ) into (IV-4) 



5 



"\ 



[c ('em)] [a a] 



r If ArOnd)+ 5 



(c^cm)][AA3 'V/’md) 






Ccfcm)][AA] 



DC 



(ir\d)lo^^ 



which is the mechanization equation in terms of ( DC ) ( i n d } ( l on g ) 
and A r (i nc j). Further substituting (IV-1) and (IV-2) into (IV-6), 



(IV-6) 



( 1 1 


A. ^[c(cm)J [A A) 1 


$ 1 1 

| '- , rc(crv\)]CAAl ’ 


\a / 2 S 

^ me s '[c(cm)3CA/\J 4 


'Vftru.e) ~ r»r(true) 

^Lcfcm’DrAA] 


S 5 

^tc(cm)3[AA) v,/ [c(cm)Kaa]> 



/ A r (trtte)“ 



=(C)A r -f- 



^[C(C rr.)] CA Al (t)A r 
[a a] 



-f 



XC(crYi)3[AA] 



(C)A 



5 



[cfc w>vj[aaO 



( lv-7) 

This control equation indicates that the value for S[c(cm)][Aa1 

should be unity, and it will be so chosen in further considerations 
of the problem. Even with the position error eliminated by the 
proper selection of S[(c)(CM)] [aa] * there remains an error caused 
by the rate of change of A r ( true ). Since ^ r ^true ) a me asure 

of the distance travelled over the surface of the earth, and because 
airspeed is constant, 



A 



r(tru«) 



= V. 



Ce-f'ai*-)! 



( IV-8) 



JO 



This means that the steady state error is a function of 
wind velocity, and that this equation is therefore unsatisfactory. 

B. In order to improve this response, the first integral of 
A r(app) added to the second integral in the preceding mech- 
anization equation. The resulting equation is 




dtdt +- 5 



[cicmO Caa3 




from which, with (IV-3), 



^[cfc^)][AA] ^r(\r\6) '^Cc( , cm)3[A A] ^ 



Further substituting ( IV— 1 ) and (IV-2) into (IV-10), 



i S 



Lcfc^)D[AAj 



w 



NE 



Ec(cm>3LAA] 



A r 



(true) 



- I- 



W, 



we 



[A A] 



• I 



A 



. . = (C)A r + 

r(true) C 



^[c(cm)3[AA3 



[C(cm)J [A A] 



In this control equation, there is a forced error resulting 
from the acceleration of the wind which can be removed if 
1 ? 

7 i — is made equal to W^ne* Under this condition, a 

S [C(CM)] [AA] 

stable control equation results with a forced error resulting 
from no effect of the wind of lower order than the time rate of 
change of acceleration. The homogeneous equation contains damping, 
but has a period of about 64.4 minutes, because the natural fre- 
quency of the system is fixed as W 2 NE when the value of S[ C ( CM )j [*aa] 



is chosen. Bee , use it is desirable to increase tnis natural fre- 
quency, in order to reduce the time during which transient errors 
are of importance, other mechanization equations are examined in 
this cnapter , in an effort to acquire control over period as well 
as damping in the control equation. 

C. The third mechanization equation to be considered uses 
the third integral of A r ( app j to move the controlled member, by 

operating on A r ( ind j and its first tnree integrals. This mech- 
anization can be expressed 




dtdtdt +■ 5 




r(\ nd ) 



dt dt+ S 



[c(Cm)J [AA - ] 




rbnd) 



dt + 5 



[c(: m )3 



Differentiating equation (IV- 12 ) three times, and substitut- 
ing equation (IV-3) into ( IV-12 J, 

* * 

• « r » 

^r6nd) ^ CC(C r* )] [A A] ^iVi'nd) S Cc(cm)][A a] ^r(ind) (OC) <M „f 



f'urther substituting equations (IV-1) and (IV-2) into eq 



(IV-13), 



-A 



y ( i rue) 






Cc(cm)]CA A3 



Tcfc^OCA Aj 



w :5 



A 



5 



Cctemi][AA3 



r(irae) 



me ^ [c(c^'>3Caa3 



[cfcm)] [A A3 



A 



f-trup) 



This control equation again contains a forced error propor- 
tional to tne velocity of the wind, and when the proportionality 
factor ^ [\, * j;, . ; ] £ AA ] -*- 3 made zero, tne homogeneous equation becomes 





(IV- 13 ) 



(C) A r -h S ec(c^CAA] .. ■Otc^oCAA ] (c)Ar + 



[c(citi')3 Ca a3 



5 

5 



(C)A r 



CC(cmOC A A3 



5 



CcfcmOCA a] 



( tv-14-) 






0 

unstable. 

D. In an eifort to improve this situation, the second in- 
tegral of A r ( a ppj is added to the third integral in the preced- 
ing mechanization. Then, 





dt dtdt “h 5 [cfcm)] [a a] ff dt dt + 5 






Cc(cm)J tAA} 



/Arf.nd) dt +5 



CcfcrviOtAA] 




... A 



(app) 



dtdt dt 4-5 



CcCcmOCAAJ 



A r(ap p; dt dt 



from which, with eq (IV-3), 



( iv-15) 



Ur-.- Ai . .*:* . - A 



Li'.j ^ j b.. . ~’p. (OC) 



= rock 



(i*v4) lor»0 



Further substituting eqs (IV-1) and (IV-2) into eq (IV-16), 



7 1 ^rc(c^TH:AA3 



I /t c 

W we ^ccfcmOdAA)^ 



• M 



A 



S 



'"(•true) 



_ J^racmUCAAl 



[CfchnBCAA) 



W 2 5 ... , 

NE ^CcCcmlDLA A} J 



A, fc j=fcfi r + (C)A + ^^ (C)A r + 



CC)A r 



Ccfc w)]CA Al 



Ccfcitvl] [ a Al 



Cc(cw>)DLaa1 



( IV-16) 



( iv- 1 7 ) 



The lowest order of the forced error in this control equatior 

is that of the acceleration of the wind. The coefficient of this 

terra, however is ^ Jt( . CM)I [AA] 1 1 

s CC ( CM )3 [ AA] W 2 NE s tC(CM)][ AA] J 

s fC(CM)] [AA] is chosen as equal to -1 — , the acceleration .term 

w2 NE 



disappears and the control equation indicates that the lowest 
order of forced error is that resulting from the time rate of 
change of the acceleration. The coefficient of this term is of 
similar form to that of the acceleration term, and is 



S, 




S[C(CM)] [AAl 
W *NE S [C ( CM)] [AA]. 



It thus appears that the acceleration 



rate term can be made to disappear if 
as W 2 ne . 



S[C(CM)][aa] i3 ch03en 
S [C(CM)][ AAl 



From Routh T s stability criteria for a cubic, the sign of 
each term in the homogeneous equation must be oositive, and 



c 




^[c(cw)Kaa1 


> 1 


^[CfcmDCAAl 




5 

wcmncAA] 




5 


^ 5 

^[C(Cm)]QA A] 



Equation (IV- 1£) shows that, with S [q ( qjvj)] [ aa] * — » the system 



NE 



is stable only if — A-i . i 3 greater that w2 NTr . 

S [C(CM)][AA] ne 



The 



acceleration rate terra can therefore not be eliminated if the 
system is to remain stable. The equation that has been considered 
here, if it is physically possible to mechanize it, appears to 
provide a satisfactory system. It remains to be seen whether 
additional terms in the mechanization equation can improve the 
response. 

E. Consider next, therefore, the performance which results 
when to the mechanization equation just considered, the first 
integral of A r ( a pp) is added to the second and third. This mech- 
anization can be written 




A-omi) dt dt dt + 




v(\rs<\) dt dt ^[cfc^CAA] /A-(md) dt 



CcfcrvOJ [ aa] ^rfind) 



( IV- 18 ) 



Arfapp) 



dtdtdt + S 



Cc(cm)] CA/U / / ' 't-r^pp) 



A„*.„„, dt dt “ 3 Cofc vn«CA A3 / A r( - ap p) dt ^’ V 
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H 



from which, with eq (IV-3), 



S [c(cm)](XA3 A r(inch = (DC), : +S 



(ind) lon^ tc(cwl 3 t ^ A 3 



• • * 

CDC)^ n<J ) \ on y~\~ S[c(cm)3CAA) ^C) (|nd)loo3 



(IV-20) 



Further substituting eqs (IV-1) and (IV-2) into eq (IV-20), 



1 5 



CacmQCAAl 



W l 5 

NE 



^ r (t rcoe) 1 



Cc(c^3[aa3 



Lc<cm) 3 C a A 3 



W 2 5 

* NE ° 



[c(c^3Ca a] y 



, i 

A , -\ 

^V(tru.e) 



5 



^“fC)A f +^“(c)A r + 



ME ^cc(cm33 l a aJ 



5 



rc(cm)3CA aQ 



5 



Cc(cmOC a A3 



Cc(cm)] La’ a] 



(C) A, 



( IV-21 ) 



Here again, the lowest order forced-error terra contains the 
acceleration of the wind* The coefficient of this term is now 

■- j~ - , and does not contain any controllable sensitivity* This 
W NE 

ayatem is therefore less satisfactory than that of par 2-D* 

F. For the sake of completeness, A r ( a pp) is next added to 
its integrals in the mechanization equation just discussed. This 
gives 




A, W) dtdtdt+S 



[C(cm)}[A.A] / / ^r(md) ^ ^ ' ^[c^m)DEAA 3 




r(A pp) 



dt dt dt JJA 



r(App) 



dt dt ~\~S r , dt H~ S Cc(cm):i [A A r(app ) 



[C (Cm)] LA A] / r ' , r(z. pp) 



( IV- 22) 



from which, with eq (IV-3)> 

. •• 

^t(cw)3CAA3 ^^ihd)loA^^[c(c^3CAA3^^'6hd)\on3^ 5 



LC (C Yvrtl L A A 3 ^ r^a ) lo nj 0 
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Further substituting eqs (IV-1) and (IV-2) into eq (IV— 23 ) 



W 



A -f 

'Vftrue) 



«• ** 



I 5 



ME 



w, 



*L 

WE 



Cc(cmQ[AAl A L 

r(W W15 



Cc^cmnCAAl 



Ay- (true) 



Ccfcm)}LAA] 



NE ^Ccrcm)3CAA3 



W15 



A. (w -fcfi r + ^»"' ;ci -Y C)A r + 7 £<g !iaa (C)A, + 



NE ^CcfcrnUC A A] 



5 



Cc(c m)] [ A A] 



Cc(cm)3C a’a] 



5 



fc)A r . 



CcrcmnC’AA] 



This system, like that of fa r 2-E, contains an unremovable 
steady-state- error proportional to the magnitude of the wind 
acceleration* 

G. When examining possible mechanization equations, it 
next seems reasonable to add to the system of par 2-D the fourth 
integral of A r ( i n d ) » usin S the fourth integral of A r ( ap p). This 
gives 



A,. dt dt dt dt + S 



from which, with (IV-3), 



[cfcmVftAA] 



Arfind) dtdt ^t+'S [cfcm q C j; A] / / A rfin j) 



[cfc mV] C A A] 



dt dt + 5 



Ccfcm)][ 



A A] dt “^^[cfcmULA'A] ^r(ind) 



A r ^ pp) dt dt dt dt 



• ' |4 I \ 

^[c(cmv][AAj A ^Cc(cml3CAA3 Ar(inil)^* ^[c(cm)]CAA] ^ccfcmfllA A"] r(i'nd) ' '(indl lohj 



Further substituting eqs (IV-1) and (IV-2) into (IV-26), 



♦ 



( IV-24) 



(IV-25) 



( IV-26) 
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^r(' true) 



ic(cm)} C a A3 



^[c(cm)] C A- 



WWW 

A — < 

'V ft: rue) 



A] 



Zc(c rn)3rA/v) 



I 



Ccrcm^rAA] 



Wl 5 



HE ^[CfcwflCXA]' 



U ~ 5 - 

(true) o 



■ £=atiil \ M =(C) ('C)A,+^ ffi = atii (C)A r + Cc)A, + - 

^CcCcm^L A*A] "^CcrCttiVJt'X'AJ ^Cc(cm)3t*A A3 ^ 



Cc(cm)3C A A] 



Cc(crY03C’AAj 



(OK 



(IV-27) 



As was true with the corresponding system of par 2-C, an 
error caused by the velocity of the wind results. 

H. When the third integral of A r ( a pp) is added to the 
fourth integral used in par 2-G, the equation becomes 



^jj J ^ ^ dt + ^[ic(cm)3CAA] ind) dtdt dt 4~ S Ccfcm u C A A 3 dt dt ^5 CC (cmflcXA] 4" 5 Cc ( cm 3KAAl A r (/ h d)' 



CcfcnnGC A A3 



A K4pp) dt dt d: dt + S [c(cmnc/U) /// A r(lpp , dt dt dt 

(IV-2b) 



from which, with eq (IV-3) 



^Cc(cmTJCAA3 ^r(fhd) - ^ ^LC(cvv03[*a A3 ^rfind) ^ ^Cc{cm)]C*A a! ^Kmd) long'd"" ^c^(cm)3CAA3 ^^0(ihd ) long 



(IV-29) 



Further substituting eq (IV-1) into eq (IV-29), 



• • »• 



A J Cc^cHaLXAl 

rf-true) J Q 

^rc(C»Y\XIt*A'A3 



Cc(cw)3:aa3 



W15 



^ (true) 



HE ^CcfcmtftT/yJ. 



TCfcmnCAAl 



CcfcmnrX'A] 



W" 5 

Vy HE 



Cc^cm)]t’A a3j 



A. (M -(C)A,+^ ti£ “ fc)A r + -|^aa(C)A+ %***(C)k r + t 

O r -/ w r ./ Jr- 



(ClA, 



CcCcrnHC A A3 



Ccfcm)]C a /*\) 



CcfcmflC X A] 



Ccfcm; ] c A a3 



(iv-30) 



Here, as with the system of par 2-D, the control equation 



indicates that the lowest order of error in the steady state 
results from the acceleration of the wind. The coefficient of 



thiatarmis S CC(CM,)][ A>3, _ 



S[C(CM)] [Ta] w2 NE S[c(CM)]Ua1 
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k 



< »; 



'* 11 / 




As before, if S [c(CM)l [AA] is made equal to — X— , the 

w NE 

acceleration term disappears. The coefficient of the acceleration 

rate term is _ !CC(CM)]LAA: ^ . It appear3 

s [C(CM)]lTa] w2 NE s [C(CM)]i AAl 

that this term can be made to disappear if _ LC ( CM )] C AAl i s 

S [C(CM)] [ AAl 

chosen equal to — JL - 

W^NE 

From South's stability criteria for a quartic, the sign of 
each term in the homogeneous equation must be positive, and 



S ... S .. S . _ K 

v *^CC(cn' 03C AAl A} LC(CmU[AAj v r ^,„nrT*i W, 



Ccfcm)] C' a'aJ |^Ccfcin]CAA 3 



S ) > 0 

^[cfcrrOJtAA] ' 



( IV-31 ) 



If S ^C(CM)][AA] is made ec l ual t0 - o 1 - > ec l (IV-31) becomes 



W 2 



NE 



(r ^ 




[ (V 


J^[c(cm)] [A A] 


.J 


1 ^CcrcmDCAAj 


5 

yCcfcmOCAA] J 


II 


r 

£ 

z. 

m 

S , 



V 



^ccrcmd C’a ad \ ^CcfcmUCAA} j 



s 



[cfcm)]C A A] 



UC(c C A Aj j 



1 r 



< 



5 



[cfcm)JCAA] 



>0 



(IV-32) 



Therefore, in order to eliminate the acceleration rate terra, 

- S [C ( CD )] [*A*Al [ s [C(CD)Haa]| 2 raust be greater than zero, which 
is possible only if S ( Qp j] [*^ 4 ] becomes negative. Routh's first 
criterion of stability, which forbids variation in the signs of the 
quartic coefficient, is then violated, and the system is unstable. 



vr 



ONFiDE 



' f 
/ i 



f . • 

* i( 



\ 



Once again, therefore, it has not been possible to eliminate, 

In tne steady state, any effect of wind of higher order than 
acceleration. Mechanization eq (IV-23) appears to give the 
same accuracy, in the steady state, that results from the use 
of eq ( IV- 15 ) • In choosing between these equations, it is 
necessary to compare dynamic response and ease of mechanization. 
Before an examination is made of this problem, a study will be 

made of other mechanization equations. 

4 

I. In an effort to eliminate the effect of acceleration 
rate, the second integral of A r ( app ) is next added to the mech- 
anization equation of par 2-G. This gives 



Ijfh r(.nJ) dt dt dt dt + 5 tefc(B)3ti JJJ, A r(md) dt dt dt +5 



Cc(c*m)3C a*: 



dt dt +5 



IcUmntfftj ‘ \c>rA) dt dt A r(in< j) ^/^^A rfap p-, dt dt C(t dt +5 



A,,, dtdt dt -H 5, 



from which, with eq (IV-3), 



CcfcmHUA]// / / ~Vipp) Ul T~ J K fcml 3 CAM , 1 A r(tff) dtdt 



(iv-33) 




Cc(cm)3CAA3 ^ Cc(cm)] LAAO (tnd)lon^ ^ £c (em’OC aaD lon<j ^ CcfcmVj A3 



• » 



further substituting eqs (IV-1) and (IV-2) into eq (IV-34) 



_ J I ^CctcmPCAAl 
1 1 v r 

^CctewOCA’AI 



A / — < 

r(tru.e) ' 



'cc^Kaa] 



tc(cm33 CAA] 



CC(cmnC’A*A3 



W 1 5 

VV NE '"'CCfcrwTKA A] ^ 



Arftrut)’ 4 " 



I 



^WE S £c (c A* Al 



A,i,^=fc)A r + |“ Mlll fc)A, (C)A,+ ^ ,1 - B1 “ 1 

Orw^nrTo ^CCtCfnUCVXl ^ r - 



(c)A r H- 



CC(crrv\3CA A] 



:c( C r«m A A] 



5 



fC)A t 



[cfcm)K*A A3 



(IV-34) 



(iv-35) 



This control equation shows that, instead of improving the 
response, the use of eq (IV-33) has resulted in an acceleration 
error which it is no longer possible to eliminate. 
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J. A final effort to remove the acceleration rate terra is 
made by adding the fifth integral of A r (*pp) and of ^r(ind) to 
the mechanization equation used in par 2-G. This gives, 




A r(Vnd) dt dt dt dtdt-5 



from which, with eq (IV-3), 



(C(C ">)] [ a A] 




A Ki ^, dt dt dt dt -hS 



Ccfc ">)3£ A A) 



'A 



ir»d) 



dt dt dt + 5 refcw03ffA3 //A^ md) dt dt+5 



A At + S A 

CCCOv)jf AA^/ U ”a”a 3 w'l^d) 



'A r , Jpp) dt dt dt dt dt +-S CeftfB0[AA , 



A r( s.„, dtdtdtdt 



(iv-36) 






CUo^XH A a] ^CcCct^) 3 C A A"] ^ Cc (c m)} C A A 3 hd) 



(lv-37) 



Further substituting eqs (IV-1) and (IV-2) into eq ( IV— 3 7 ) » 



— A — 

'' t (trur) 



CcfcmflC a’a 3 



A 



Lc(c»^)K a A 3 



L c(c C A A 3 



C ' *rftrut) 

^cefc^nc a a3 



5 



Cc^cw) 3 C’a'a 3 



\J x S 

^ NE ^CcfCrnUC’ a‘*A 3 



A — 

r Wtrue; 



Cc.(cm^3CA K] 



Ccfc "M 3 C'a"a 3 



I 



W l $ 

we '■'oefc^nc 



-k w - £) A, + 1 “=“^ COA,+| 

a Ajy 



CCf'cm^C A A 3 



fc)A r + / c(wl3[ '^ 1 (t)A r +r^^ (C)A r + 



Cc<'cnM 3 ["i’' a 3 



Ccfc^nt. a a 3 



5 



Wc»))t A A] 






Cc^^OC'a A3 



5 



(C)A, 



[c l A a] 



1 1 v-38) 



The coefficient of the wind acceleration term of the forced 



error now has the form 



S CC ( CM )] C AA ] 

s Ic(cm)][Ta] w ‘ne S tC(CM)][ _ A A] . 



7 



M|*4 



To eliminate this, 5 [c (CM)][ *aV) roust be made equal to — ^ — » 

W NE 

before. In this equation, the wind acceleration term has the 



as 






coefficient 



S [C ( CM )) [ AA] _ s tc ( CM )3 C AA ] 

S [C(CM)][Ta] w2 NE s CC(CM)][Ta] 



, If this is 



to be eliminated, 



S EC(CM)][ AA ] 
L S tC(CM)l [AA] 



must be made equal to — ^ 



W‘ 



NE 



< 
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r r r. ( 

L *- j •-} 
< 1» i ^ 



* i / . 



Using Routh's stability criteria for a quintic, the co- 
efficients of the homogeneous equation must all be positive, 
and 



irrcm)]riA] 



[cfcmflt "a*a] 



V . ^ t «»• 

Cc(cm1][ a A] 




>o 



The system is stable, with proper choice of the other co- 
efficients*^ S [C(CM)][ AA] * s mac ^ e equal to — , eliminating 

w NE 

the acceleration term* If this substitution is made in eq (IV-39) 

1 



s ••• 

and if in addition the substitution — [C(CM)][ AA] 

S [C(CM)H AA] 

made, 



W‘ 



is 



NE 



(IV-39) 



Tc(onrj[AA] 



5 



rcfcrvpjt a a3 



Cc('cm)3C X A] 




> 0 



( IV-40) 



This equation simplifies to 
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5 



Cc ('em!] Caa 3 



5 



LcfcrrODt A aJ 



5 



Cc(cir03 C 'a A3 




S S . -S 

"CdcnflrAA] "^Cefem)] Ca a] [c fc ml] C a a] 



'-\cfcmiK a" A] 



- s 



CcCcm)3CA a] 




( IV-41) 



Since all of the sensitivities are positive, this equation 

2 

teW(CM)] [7a] 

_ •••• 

S [C(CM)][ AA] 



can be divided by 



, provided that 



only S Jq ( Qp y\ j is maintained smaller than ( ■ 

\ ^ wE 

without changing the inequality sign. This yields 



5 



C. — c. 

Ccfcm^tAA] ^[c(cm)3LAA] u 



[CfcrrOJS; A A] 




C*A*A3 






HJ 



5 S 



ccfc^ni: a'a3 ^{.cic m3 Ca a] ^[c(fw)][A a] 



-s, 



or, simplifying, 



-S 



rc ('em)] C A A3 



> 0 



Since ( CM )] [**A A ] must be kept positive in order to satisfy 

the first Routh criterion, both criteria cannot simultaneously 
be satisfied, and the system becomes unstable if an attempt is 
made to eliminate the acceleration rate steady-state error. 



( IV-42) 



(IV-43) 











t 



3 . Implementation of Equation s 

A. Four mechanization equations have now been discovered 
that eliminate steady-state error caused by the acceleration 
of the wind. It has not been possible to eliminate any higher 
order terms, but is seems unlikely that any steady-state error 
resulting from the acceleration rate of the wind can have appre- 
ciable effect. The simplest of these equations was eliminated 
because natural frequency could not be controlled. The three 

i 

usable equations are (IV-15), (IV-22), and ( I V-3 6 ) . In eq 
(IV-15), two of the three sensitivities can be controlled; in 
eq (IV-28), three of four sensitivities can be controlled; and 
in eq (IV-36), four of five sensitivities can be controlled. In 
order to obtain the most satisfactory possible dynamic response, 
the use of one of the more complicated equations may be indicated, 
but this choice of equations is subordinate to the considerations 
of ease and accuracy of mechanization. A practicable method of 
implementing eqs (IV-15) and ( I V-2 B ) is now presented, and the 
complexities compared. 

B. In mechanizing eq (IV-15), (DC ) ( i n d } ( long ) is usec * 
instead of A r ( ap p) as the output of the pendulous accelerometer, 

as explained in par 1 of this chapter. The resulting eq (IV-16) 
is restated here for convenience, in integrated form: 
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( IV-44) 
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Now, an imaginary auxiliary direction is established, as 
shown in Fig IV-1. This direction is controlled, with respect 
to the reference vertical, by the equation 
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The controlled member is positioned from the reference vertical 
through the use of the auxiliary direction by the equation 



(IV-45) 



A Kind) S Ccfc~)fauxV) [aa] / / (ihcOlow^ dt dt-S Cc(C>*VdwX)3[ a A 3 j ^‘r( i U -**.) 



( IV-46) 



These equations, when combined with the relation 
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This is identical with the desired eq (IV-44), with the 
relationships among sensitivities that 
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FIG.1Y-1 

AUXILIARY DIRECTION FOR CUBIC MECHANIZATION 

EQUATION 
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Equations (IV-45) and ( IV— if 6 ) can be mechanized easily, as 
shown in Fig IV-2. The mechanization of eq ( IV —15) is therefore 
feasible. 

C. Equation (IV-28) is now to be mechanized in a manner 
similar to that used for eq (IV-15). The integrated form of the 
equation, using ( DC ) ( i n d } ( l on g ) instead of A r ( ap p), as before 
becomes 




For this equation, a pair of auxiliary directions is esta- 
blished, as shown in Fig IV-3. These directions are controlled, 
with respect to the reference vertical, by the equations 





dt dt 
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FIG. IV-3. AUXILIARY DIRECTIONS FOR QUARTIC 
OR QUINTIC MECHANIZATION EQUATION. 
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The controlled member is positioned from the reference vertical 
through the use of the auxiliary directions by the equation 
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These equations, when combined with the relations 
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This is identical with the desired eq (IV-52) with the 
relationships among the sensitivities that 
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Equations (IV-53), (IV-54), and (IV-55) can be mechanized 
almost as easily as were eqs (IV-45) and (IV-46), as shown in 
Fig IV-4. The choice between the cubic and quartic mechanization 
equations is largely, therefore, to be determined by a comparison 

of the dynamic responses of the system when these equations are 
used* 

D. The method of mechanizing the quintic equation, ( IV— 3 6 ) , 
13 similar to that for the cubic and quartic equations just dis- 
cussed. In integrated form, using (DC) ( ind ) (long) instead of 
A r(app)i this equation becomes, 
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MECHANIZATION EQUATION. 
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mechanization equation. 
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Tne two auxiliar. directions established by eos ( I V— 5 3 > 54 , 5 b,and 
57 ) are required to implement this equation. 

The controlled member is positioned by the use of these 
equations and the equation 
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When these are combined, the control equation becomes 

\ 



( IV-64) 
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This is identical with the desired eq (IV-63). The mech- 
anization of the equations from which ( IV— 65 ) is derived is 
shown in Fig IV-5* This mechanization is, of course, more 
complex than the systems previously discussed. Although the 
mechanization is entirely feasible, the quintic equation will 
not be considered further unless the dynamic responses of the 
system using the cubic and quartic control equations show 
themselves to be unsatisfactory. 

4. The Closed-loop Syste m 

In order to examine the two mechanization equations between 
which choice is still to be made, it is necessary to examine the 
response of the closed-loop longitudinal system to wind forcing 
functions, using each of these mechanization equations. Response 
equations for particular mechanizations have been derived. The 
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general closed-loop equations are derived in this section. 

The general mechanization equation can be written, for 
this study, as 
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The linear distance of the missile measured along the track, 
from the initial reference point, can be written as 
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These relationships! with eq3 (IV-1, 2 and 3) and the small 
angle pendulum equation 
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combine to form the closed-loop diagram of Fig IV-6. 
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FIG. EE-6. RANGE INDICATION SYSTEM CLOSED-LOOP DIAGRAM. 
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The initial conditions in eq (IV-67) do not affect the 

e 

response of (C)A r to V Qg _ ( a ir)] > an< * can be t & ken aa zero. 
Figure IV-6 can then be simplified to Fig IV-7. 

In accordance with servomechanism theory, the overall 
performance equation of the range indication system can be 
written 
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mechanization equations to be examined can be written 
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(IV- 72 ) 



Using these equations, numerical values can be assigned to 
the sensitivities, eq ( IV-69) solved for various types of wind 
forcing functions, and the responses using eqs (IV-7 2) and 

(IV-71) compared. Numerical values are first assigned to the 
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FIG.nr-7 SIMPLIFIED RANGE CLOSED LOOP DIAGRAM. 
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cubic mechanization eq (TV-7\l), and then to the quartlc eq (IV ”72). 

5. Selection of Numerical Values 

A. The homogeneous control equation for the cubic system can be fac” 
tored into a first order and a quadratic term, becoming 

(f> + — )(p 1 +2(DR)Wp+W I )= 0 , (iv -73) 



which expands into 
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Comparing coefficients term by term, with 
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Then, for any desired characteristic time and quadratic natural fre- 
quency, a damping ratio and set of sensitivities can be computed. Notice 
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that, once the characteristic time and the natural frequency have been 
selected, no control remains over the damping. 

B. The homogeneous control equation for the quartic system can be fac- 
tored into two quadratic terms. As in the track control system, these are 
made equal, giving, 

(p 2 + 2 (DR)Wp +W 2 )(p l +2(DR)Wp +W> O (IV-78) 



This expands into 
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Equating coefficients with those of the quartic control equation, with 
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Plots follow for responses of the cubic and quartic control equations, 



with various sensitivities. 
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TIME IN MINUTES 

fig. nr- a. response of cubic range indication system to 94.2 mph. rotating 

WIND OF 30 MINUTE PERIOD. 
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FIG. nr- 9. RESPONSE OF QUARTIC RANGE INDICATION SYSTEM TO 94.2 MPH. ROTATING 

WIND OF 30MINUTE PERIOD. 
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TIME IN MINUTES. 

FIG. I2-I0.COMPARISON OF CUBIC AND QUARTIC RANGE INDICATION 

SYSTEM RESPONSE (DR) = 0.295. 
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CHAPTER V 

CONCLUSIONS AND SUGGESTIONS 

1. Disturbances to the Control Systems 

In examining the success of the mechanization systems whose responses 
have been analyzed in the preceding chapters, it is necessary to consider the 
types of disturbing inputs which will affect the systems. These disturbances, 
or winds, are of three major types. First, there are gust disturbances of 
fairly high frequency and short duration. These are oflormously attenuated 
by the long time constants of the systems analyzed, and should cause no 
trouble. Second, there are acceleration pulse type Inputs, such as would be 
encountered when crossing weather fronts, where the wind velocity changes 
very rapidly to a finite and fairly constant value. These can be analyzed by 
the step function responses of Chapters m and IV. In Fig V"l, the path 
through the "eye" of a hurricane has been approximated with step functions 
of acceleration, and in Fig V”2, the response of the track control system is 
shown to be satisfactory, even with this violent disturbance. 

There is, however, a third type of disturbance; that in which a wind of 
approximately constant velocity slowly rotates.* This condition has been 
approximated by sinusoidal inputs as shown in Chapters in and IV. The error 
if the sinusoidal forcing function continues for a long period of time, is seen 
to be very large. A wind of about forty knots, rotating with a period of about 
one hour, will introduce, in the steady state, deviations from the track, in 

the track control System, of about 20 miles. 

✓ 

•Ref. A. G. Bogostan Notebook dated March 18, 1947. 
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FIG. V-l. APPROXIMATION OF A HURRICANE USING 

ACCELERATION PULSES. 
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TIME IN MINUTES. 

FIG. -I- 2 RESPONSE OF TRACK CONTROL SYSTEM TO IDEALIZED HURRICANE. 
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Three ways of correcting this condition suggest themselves: 

a. Build systems with variable coefficients, which provide the neces- 
sary high damping with step-inputs, and have a different response to such 
low-acceleration inputs as sinusoidally varying winds. 

b. Increase the order of the track control system so that the forced 
error results only from high order greatly attenuated terms. 

c. Modify the system to reduce the effects of the forced errors. 

These three methods are all examined in the next section of this chapter. 

2. Suggestions for Improvement of System Response 
A. Track Control 

a. A variable coefficient system can be of value only if changing the 
damping and period of the system will improve the response. Since the co- 
efficient of the acceleration rate forced error term increases rapidly with 
increase in damping, it appears that this should be possible. A series of 
track control responses using the same magnitude of disturbance but dif- 
ferent damping ratios has therefore been examined, Fig. V-3. It is seen 
that no choice of damping is satisfactory. If damping is entirely removed, 
the forced error of the simplified system reduces to zero. Actually, a 
forced error will still exist if there is any vertical wind acceleration, or if 
non-ideal components are used. Then, if the frequency of the rotating wind 
is equal to that of the system, very large errors can result. This approach, 
therefore, does not provide the needed improvement of response. 

b. The purpose of this section is to investigate the response of a 
quintic performance equation. From eq m _ 26 the performance equation 
may be written: 
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For this study it is assumed that the denominator of the performance 
equation may be factored into the form: 
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Then expanding and equating coefficients of like items, 
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Examination of eqs (V-2, 3, 4, 5, and 6) reveals that if K and (DR) are 
selected, then CT and all the sensitivities are fixed. General expressions 
are now developed for determining the numerical values of CT aadul certain 
ratios of sensitivities. 

From simultaneous solution of eqs (V-3 and 5) 
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From eq (V~6) 
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From eqs (V - 5 and 7), 
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From eqs (V-4, 6 and 7), 
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From eqs (V“2, 4 and 6) 
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Limited time permits the evaluation of only one numerical solution of 
the performance equation. 

Let K = 2 and(I®)= 1. 

Then, 
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Substituting these values into the performance equation and using as an 
input Vj-g _ (air)] = ft/min^ sin t, 
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(C) A (Vert) (tc) = .8 minutes of arc (.8 nautical miles) 

Comparing this answer to those obtained in the plot of results in Chap- 
ter in indicates that the quintic performance equation produces a better re- 
sponse than the quartic one. It should be mentioned that no effort has been 
made to select optimum parameters for the quintic. 

c. It is seen, by an examination of the track control system perform- 
ance function, that no control over the gain setting exists. As was shown in 
Chapter in, it is desirable to have the sensitivities negative for the quartic 
performance equation. Making this change, eq in-25 becomes; 
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where sensitivities are now positive numbers. 

If a controllable gain coefficient could be introduced into the numerator, 
so that an attenuation of the forced error could be provided, then the system 
might be made to perform without excessive errors. This could be done, it * 

appears, if the wind acceleration term could be attenuated after entering the 
system. This cannot be done, but an approximation of this desired result 
is achieved if additional feedback loops are provided as indicated in figure 
V-4. The performance equation for this system, when (PF) ( aer0 ) (i n d) is 
assumed to be exactly equal to (PF*) ( aero ) ] 
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FIG. V-4. MODIFIED TRACK CONTROL SYSTEM. 
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when the derivative auxiliary feedback sensitivity is set at zero, this becomes, 
with proportional plus first and second integral mechanization equation, 
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In this equation, a certain amount of gain control has been introduced 
Into the numerator, at the expense of introducing a term giving a forced 
error caused by wind acceleration. Houth's stability criteria show that this 
system, with proper selection of sensitivities, is stable. Limitations of 
time did not permit the authors to make a complete analysis of this control 
system. However, a sample result, selected at random is calculated. The 
denominator is first factored into two unequal quadratics, and the coefficients 



squared. They become: 
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The values chosen were 



S(aux)[AA] 

w i 

W 2 
(DR) 2 



= 0.1 
= 1.5W 

= W NE 
= 0.7 



NE 






These gave, using the equations given above, (DR), = 9.18, 

CC(CO) 1 [ A A 3 \ l~> CC(CD)KAAD \_ _- Q . .( 9[R 

=■ 0.124, (— ) = 22.84 and 



\- > [C(CD)]CAA) 



CCtCDllt A A] 



^rctco/KAAi^aoiO-M) 



20.72 



The magnitude of the maximum forced error was found to be 28.5 nautical 
miles, and the disturbance which gave this maximum deviation was a sinu- 
soidally varying wind with a velocity magnitude of 265.5 nautical miles per 
hour. This is a considerable improvement over the uncompensated system. 
Next, was set at unity to remove the accleration term, and 

S (aux)[AA ] was set > at a ranc * om va * ue * 

The performance equation, with quartic mechanization, becomes 
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Equating coefficients, as before, to those of a quartic with the desired 
characteristics, and selecting, at random, the following values: 



*(aux)[AA] 



= 20 



W 1 - W 2 = 2W NE 



Answers were obtained making 



CC(C0)][A 



cccconcA 



S A]/ 



(DR) 1 =(DR) 2 =2.75 > 1 CcCCP)]CAA] | s 4.78 and) 

S>CC(CD)][AA] 



= .0885 

9ir 



, [CCCO)][AA] V C(alr)-M] 



-3.25 



With the same wind input as that used in the preceding example, the maxi- 
mum error was found to be 48.8 miles. This, again, is an improvement 
over the uncompensated system. 

It appears, since in the first example the error was caused chiefly by 
the acceleration term, and in the second example by the second time rate 
of change of acceleration, that a combination of the two systems, varying 
both s ( aux )[ AA ] and s ( aux )[AA ] might S ive a satisfactory solution. Time 
prevents such a study here. A simulator would greatly facilitate the work 
of selecting optimum sensitivities. 

B. Range Indication 

Lack of time prevented an elaborate examination of the response of 
the range indication system to large ranges of inputs. Since the track con- 
trol system indicates that the critical response is to sinusoidal inputs, a 
series of solutions was made for cubic and quartic responses to a sinu- 
soidal acceleration input with a period of thirty minutes. 

As indicated in the plots of Chapter IV, the response is improved as the 
damping ratio of the quadratic in the cubic equation, or the two equal quad- 
ratics in the quartic equation, is decreased. At a damping ratio near 0.3, 
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3. TRACK CONTROL MECHAN ISM SUMMARY 



MECHANIZATION EQOAT'ON 



« 
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THIS SYSTEM HAS NO DAMPING AND BECOMES UNSTABLE IF ATTEMPT IS MADE 
TO ELIMINATE FORCED ERROR INTRODUCED BY ACCELERATION OF THE WIND , 
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4. Range Indication Mechanization Summary 

Examining the control equation obtained in Chapter IV, it is seen that 
it is possible to eliminate forcing function terms of orders up to and in- 
cluding wind acceleration terms. In a stable system, higher order terms 
cannot be removed. A general performance function including all mecha- 
nization equations that eliminate acceleration and lower order terms can 
be written 



(PF), 



(mech) (A A)(lon<j) 



S(n-1)T P * 5(n)T 

p n + 5 (i)t p n •••• + S (n _ z)T p 2 (V_19) 



If more terms are included in the denominator, velocity and position 
errors result; if more terms are included in the numerator, control is lost 
over low order terms, while higher order terms are eliminated. If terms 
are omitted, instability results. 

A block diagram showing a method of mechanizing the performance 
function eq (v~19) follows, fig (V-5). 

For both track control and range indication systems,, the response for 
sinusoidal forcing functions of 84.4 minute period appears to be the same 
regardless of the choice of damping or order of the equation in the T ' 
mechanization system. 
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TO OBTAIN THE SENSITIVITIES OF THE EQUATION WRITTEN ABOVE, EQUATE 
WITH PROPER TERM OF BLOCK DIAGRAM EXAMPLE! (for ^uaaraf »c ) 
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APPENDIX A 

DISCUSSION OF SIMPLIFYING ASSUMPTIONS 

Several simplifying assumptions were made in reaching the mechani- 
zation equations and performance equations discussed and analyzed in this 
thesis. These will now be examined. 

1. Omission of the Geodesic Acceleration Term 

The geodesic acceleration arises from the fact that navigation is 
performed on the surface of the earth, where the gravity equipotential 
forms the geoid. This geoid may be considered to be an ellipsoid with 
roughness. The deflection of the vertical due to this roughness never ex- 
ceeds seventy seconds of arc, and (especially over land masses) is usually 
far smaller. 

While the missile travels over this ellipsoid, computation assumes 
that the surface is a sphere. The shortest distance between two points on 
this sphere is a great circle. The gravity verticals fall in this great circle 
plane. The great circle plane transfers to the surface of the ellipsoid as 
a plane section, but the gravity verticals no longer lie in this plane, but in 
an arc intersecting it at the points of departure and destination. The mis- 
sile, traveling along this curved path, experiences a horizontal accelera- 
tion. The radius of curvature of this path, which is a minimum at 45 degrees 
north or south latitude, is never less than one million miles. The geodesic 
acceleration is therefore of negligible magnitude for missiles with veloci- 
ties of a few thousand miles an hour. 
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Two figures (Fig A-l and 2) prepared by Dr. W. Wrigley of the 
M.I.T. Instrumentation Laboratory, illustrate this geodesic acceleration. 

2. Compensation for Coriolis Acceleration 

The acceleration of the Coriolis arises from the fact that the mis- 
sile is moving in earth space, which possesses an angular velocity with 
respect to inertial space. This is demonstrated in the accompanying fig- 
ure (Fig A-3). Compensation is independent of the heading of the missile, 
depending only on a knowledge of the earth's angular velocity with respect 
to inertial space, which is accurately known, on latitude and on ground 
speed. An indication of latitude and of ground speed are available within 
the missile, the accuracy of which are dependent on the success of the con- 
trol system. Errors resulting from inaccuracies in these measurements 
will be negligible in any practicable system. 

A certain amount of cross -coupling of the track control and range 
systems will occur through the Coriolis Computer, but it will remain small 
if the component of missile velocity perpendicular to the track is not large. 

3. Small Angle Pendulum Assumption 

For simplicity of calculation, a linearization assumption has been 
made in the output of the pendulum. The pendulous accelerometer actually 
solves the equation 



[E - (cm) ](hor) = ^ (DC ) 

6m 



The control equations are based on the assumption that 



[E - (cm)j(hor) = 

g m 



(true) 
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FIG.A-I 

BASIC GEOMETRIC FACTORS ASSOCIATED WITH GEODESIC ACCELERATION 
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FIG. A-2.. 

RELATIONSHIPS BETWEEN THE GEOID , REPRESENTATIVE ELLIPSOID 
AND THE NAVIGATIONAL SPHERE; DEFLECTION OF THE VERTICAL. 
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. FIG. A-3 

BASIC GEOMETRIC FACTORS ASSOCIATED WITH 
THE ACCELERATION OF THE CORIOLIS 
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The validity of this assumption has been tested, by comparing the out- 
put, for the track control system, to a pulse of wind acceleration of 2700 



fied equation. The results are shown in the accompanying figure (Fig A-4). 

If necessary, a tangent pick-off can be used on the accelerometer in- 
stead of the linear pick-off postulated in this discussion. 

4. Neglect of Vertical Acceleration 



If the missile does not maintain absolutely constant altitude (or if 



the magnitude of the gravity force varies), the pendulum angle will not be 
a function of horizontal acceleration alone. This effect, together with the 
effect of an aerodynamic lag and a pendulum lag, is considered in the fol- 
lowing treatment of the range indication problem. There are changes in 
some of the geometric equations, as indicated: 

Aerodynamic lag term: 



o 

feet/min , lasting for one minute, using the transcendental and the simpli- 




(A-l) 



Accelerometer lag term: 




(CT) (pu) p +1 



1 



(A-2) 



I* (true) 




1^7 



(A -3) 
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FIG. A-4 EFFECT OF LINEARIZING PENDULUM OUTPUT 
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FIG. A-5. INTRODUCTION OF PENDULUM LAG 

IN RANGE SYSTEM. 
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The mechanization chosen was the cubic mechanization equation discussed 
in Chapter IV. 

The control equation for the range indication system, using the con- 
cepts outlined above, becomes, 



(9 IR * - (cm>) (Vert)) 
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^LE'«n 
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[p 3+ Sj; c{cm3 - 1 [AA 1 P + -2cc(cm)][A A1 p + S £ C (cm)H A AT H ^ 1 R* ^[E> (cm)XVert)] l(CT)(p U ) p 4 l1 [(CT) (aercy) p * l] 



If the sensitivities are chosen so that the cubic in the denominator be- 
comes the product of a quadratic with (DR) = 0.715 and a period of forty 
minutes, and a first order term with (CT) = 14 min, eq (A-5) becomes 



a[E-(c^(Vert). T 1 (, , a [E-<cn>flCVeYfl\ , rT » a [E- (cm)](Vert^| D _ ( CT W) / _ 0 [E-(cmfl(Verm p 2 
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The effect of vertical acceleration becomes very clear. If the mis- 
sile momentarily loses lift and "drops”, the denominator goes to zero. 

The numerator, meanwhile, acquires a term which gives a forced error 
from wind acceleration. This error will change sign with change in the 
direction of vertical acceleration, but the effect of the denominator will 
keep the opposing effects from balancing. The damping in the denominator 
cubic, however, will return the missile to the track after the vertical ac- 
celerations have passed. 

Even with a downward acceleration greater than the effect of 
gravity, the system retains stability. 

Reasonable values for (CT)^ aer0 ) and (CT)^ u j will vary from a 
fraction of a second to a few seconds. These will, therefore, have little 
effect on the system. 

5. Assumption of Perfect Aerodynamic Response 

Paragraph 4 indicates that neglect of aerodynamic response terms 
can have little effect on the performance of the longitudinal system. In 
the track control system, however, with a fourth order equation or higher, 
the system appears stable if a positive pendulum angle causes either a 
left or a right rudder deflection. It seems, therefore, that even a small 
aerodynamic lag might cause the system to become unstable. The problem 
was therefore recalculated for one choice of (DR), using, instead of perfect 
aerodynamic response, a pair of equations. The first gives the desired 
heading of the missile as a function of pendulum angle; the second relates 
this desired heading to the actual heading of the missile through an aero- 
dynamic lag. 
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(A- 7) 




Z(cm>- 
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M z(cm)(ind) 

(CT)(Q e ro)P + 1 



(A -8) 



The variation of the response of this system from the idealized 
system discussed in Chapters m and IV could not be detected in the plots 
drawn by the Differential Analyzer . 

6. Assumption of Perfect Gyros and Accelerometers 

Since this system is an angle measuring system which carries as 
a reference an indication of the point of departure dependent on the accu- 
racy of the gyros, any drift in the gyros which define the vertical plane 
will give an equal angular error in range and deflection. The third gyro, 
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COORDINATE DIAGRAM AND STATIC PERFORMANCE CHARACTERISTICS FOR SI NGLE-DEGREE-OF- FREEDOM 

GYRO UNIT WITH NON -IDEAL GEOMETRY. 
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•which maintains the orientation of the track plane, will introduce a non- 
linear error. The gyros have been assumed perfect primarily for sim- 
plification of the problem. If gyros of sufficient accuracy are found to 
be impracticable, the necessary inertial reference can be found from ob- 
servation of the fixed stars. 

Permanent errors in the zero indication of the pendulums also 
introduce an error on the earth's surface equal in angle to the pendulum 
error. 

Single degree of freedom gyros and pendulous accelerometers of 
types which may prove satisfactory are being developed in the Instrumen- 
tation Laboratory at M.I.T. These units are shown diagrammatically in 
Figs (A-6) and (A-7). 

7. Omission of Constants of Integration 

Wherever possible, throughout this thesis, constants of integration 
have been assumed to be zero. For example, in the range indication prob- 
lem, the missile has been assumed to have no airspeed. This does not 
change the time response or the size of the range or track error resulting 
from wind disturbances, but will have a large effect on the distance 
travelled while the system is recovering from a large disturbance. 

At all times, the system has been assumed to be in equilibrium 
when hit with disturbing wind accelerations in the time solutions of Chap- 
ters m and IV. 
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APPENDIX C 
GLOSSARY 

The system of notation adopted for use in this thesis has been selected 
because it fulfills the following desirable objectives: 

1. It is easily learned 

2. It is adaptable to a wide range of situations 

3. It is built up almost exclusively of characters found on the 
keyboard of a standard American typewriter 

4. Any one of the compound symbols of the system is readily 
interpreted without recourse to an extensive glossary 

A few simple examples suffice to explain the operation of the notation 
system. These examples are shown in table C~3. The short table of key 
symbols given at the end of this explanation will then furnish sufficient 
information to enable any compund symbol to be correctly interpreted, and 
to provide- the necessary tools for the generation of new symbols. 

A representative list of the primary symbols is given in the following 
table. 



138 



CONFIDtN i iaL 



a 

A 

(CT) 

[(CT)PR] 

(C) 

(DR) 

(Dir) 

(DC) 

(FR) 

g 

(PF) 

R 

S 

(Sg) 
t . 

V 

(Vert) 

W 



Acceleration 

Angle 

Characteristic time 

Characteristic time - period ratio 

Correction 

Dampting ratio 

Direction 

Dynamic Correction 
Frequency ratio 
Gravity 

Performance function 

Radius 

Sensitivity 

Signal 

Time 

Velocity 

Vertical direction 
Frequency 



PRIMARY SYMBOLS 
TABLE C-l 



A 
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(aero) 


Aerodynamic 


P 


Performance operator (d/dt) 


(air) 


Air mass 


(pend) 


Pendulum 


(app) 


Apparent value 


(pa) 


Pendulous unit 


(aux) 


Auxiliary 


r 


Range 


(CD) 


Control direction 


(ref) 


Reference 


(cm) 


Controlled member 


(res) 


Resultant 


(dep) 


Departure 


t 


Time 


(dest) 


Destination 


(tc) 


Track Control; in track con- 
trol plane 


E 


Earth 


(true) 


True value 


(hor) 


Horizontal 






(Vert) 


Vertical 


(ind) 


Indicated value 


X 


Axis through nose and c.g. 


IR 


Inertial reaction 




of missile 


(kin) 


Kinematical 


Y 


Axis through right wing and 
c.g. of missile 


(long) 


Longitudinal plane; in 






range direction 


Z 


Axis through c.g. perpendicu- 
lar to XY plane, directed 


M 


Missile 




downward 


NE 


Earth natural 







MODIFYING SYMBOLS (SUBSCRIPTS) 
TABLE C-2 
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Symbol 
V [E - (cm)] 

S [C(CD)][AA] 

A (Vert)(tc) 

W NE 



A Z(cm) 

(DC) (ind)(tc) 

^■^(ind)(long) 

A r(lnd) 



Definition 



Velocity of the controlled member with 
respect to the earth 

Sensitivity for correction of the control 
direction with angle input and rate of 
change of angle output 

Correction to the angle of the vertical 
in the track control plane 

Natural earth frequency equal to 
g 

IR 

«E 



Angle about the Z axis of the controlled 
member 

Indicated dynamic correction of track 
control system 

Indicated dynamic correction of longi- 
tudinal system (i.e. range indication 
system) 

Indicated range angle 



Typical Examples Showing How Symbols are Compounded from the Elemen- 
tary Forms 



TABLE C-3 



141 



re 2 60 
A P I 3 G i 



550-DUP 
B 1 I1DERY 

C ! .MDERY 



Thesxs ^3006 

E565 Erb 

A theoretical study of 
automatic inertial navi- 
gation. 



FL 2 60 1 -*• - 1 < C[ S50-DUP 
feq 6 61 JJO CRY 






Thesis 43006 
E565 Erb 

A theoretical study of 
automatic inertial navi- 
gation. 



